IPPP/04/13 
DCPT/04/26 
15th March 2004 



Double-diffractive x meson production at the hadron colliders 



V.A. Khoze"'^ A.D. Martin", M.G. Ryskin"''' and W.J. Stirling"'^ 

Department of Physics and Institute for Particle Physics Phenomenology, 
University of Durham, DHl 3LE, UK 
^ Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg, 188300, Russia 
Department of Mathematical Sciences, University of Durham, DHl 3LE, UK 



Abstract 

The double-difFractive production of Xc and Xb mesons, with a rapidity gap on either side, is studied, 
using both the Regge formalism and the perturbative QCD approach. Due to the rather low scale, the 
exclusive double-difFractive process — > p + x + P i^ predicted to dominate, whereas the probability 
that the incoming protons dissociate is expected to be relatively small. We evaluate the corresponding x 
production cross sections at the Tevatron and LHC energies. For the double-diffractive process with proton 
dissociation, it is possible to select events with large transverse momenta transferred through the rapidity 
gaps, by measuring the transverse energy, E±, flows in the proton fragmentation regions. Then the large E±_ 
provides a scale to justify the use of perturbative QCD, and to allow a spin-parity analysis of the centrally 
produced system to be performed, by studying the azimuthal angular correlations between the directions 
of the forward and backward E±_ flows. The central production of the new X(3872) charmonium state is 
considered. 



1 Introduction 

Central exclusive double-difFractive processes have traditionally been considered as a promising way to study 
new (and old) particles in an especially clean environment, see, for example, Ref. flj. In addition they give 
information about the structure of the Pomeron and of the mechanism of Pomeron-Pomeron fusion. Here we 
are particularly interested in the double-difFractive production of, C-even, heavy quarkonium (xcXfi) at the 
Tevatron and the LHC [2131, see also Refs. Since the early days of QCD, heavy quark production has 

been a fertile testing ground of many aspects of the theory. More recently, there has been much activity in 
studying the long standing discrepancy between the NLO predictions and the data for heavy quark production 
at hadron colliders 6 . Heavy quarkonium production is a valuable tool, since it provides important information 
on the physics of bound states and, in particular, allows a test of the ideas and methods of QCD effective field 
theories, which is one of the most popular recent approaches; for reviews and references see, for example, [7]. 
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Figure 1: (a) The QCD diagram for double-difTractive exclusive production of a x(O^) meson, — > p + X + 
where the gluons of the hard subprocess gg ^ x ^^^^ colour screened by the second t-channel gluon. (b) The 
rescattering or absorptive corrections to pp P + X ^ where the shaded region represents the soft pp 
rescattering corrections, leading to the suppression factor 5^. 



Detailed experimental and theoretical studies of double-diffractive x production offer a new window on this 
topic, and may shed light on some of the unresolved issues. These processes have some surprising features. 
In order to gain insight into the problems encountered in the evaluation of these processes, we first recall an 
analogous process, which in many respects is theoretically simpler. Namely, the exclusive double-diffractive 
central production of Higgs bosons, which has been advocated as a good way to study the Higgs sector at the 
LHC 8' . Both the x a-nd Higgs production processes may be written in the form 

pp->p + (x or iJ) + p, (1) 

where the + signs are used to represent the presence of large rapidity gaps. Since the Higgs is expected to 
have mass, Mh, about 120 GeV or more, its production cross section may be calculated perturbatively, via 
the diagram shown in Fig. 1. The hard subprocess gg ^ H is initiated by gluon-gluon fusion, and the second 
t-channel gluon is needed to screen the colour. The crucial observation is that the integration over Q±, going 
round the gluon loop, is made infrared safe by Sudakov factors which ensure that the subprocess gluons remain 
untouched in their evolution up to the hard scale M///2. Provided forward proton taggers are installed, the 
process has the advantage that the Higgs mass may be measured in two independent ways: first from its decay 
products and second by observing the forward protons which accurately determine the 'missing mass' Mh- The 
predicted event rate and the signal-to-background ratio make such a Higgs search feasible at the LHC, but 
unfortunately the rate is estimated to be too small for this process to be of use at the Tevatron . 

On the other hand the rate of double-diffractive x production at the Tevatron is expected to be much larger. 
Indeed preliminary studies indicate that such events may have already been detected ^HI- Unfortunately such 
data cannot provide a reliable check on the perturbative prediction of double-diffractive Higgs production. 
Unlike the Higgs, the Xc and to a lesser extent the Xb) are not sufficiently massive to make the double- 
diffractive cross section predictions infrared stable. Moreover, the predictions depend on the gluon distribution 
to the fourth power in the region x ^ 10^^ and ^ 1 GeV^, where it is not sufficiently reliably known. 
Nevertheless, motivated by the forthcoming experimental data, and by special features of the production process, 
it is informative to study double-diffractive Xc and Xb production in more detail. In addition, the interest in 
double-diffractive Xb (and even Xc) production has recently increased, since these processes may be measured 
at the LHC with luminosities of 10^" - 10^^ cm-^sec"^ by the CMS and TOTEM collaborations, with the 
special optics of the TOTEM detector m. These forthcoming, and planned, measurements will sharpen our 
understanding of the physics of these types of processes, as well as helping to find better ways to select double- 
diffractive events. 

Since M-^ is not large enough to justify the use of perturbative QCD, we will start our discussion using a 
Regge framework. We show how this formalism can be matched to the perturbative approach, and we will use 
perturbative QCD to evaluate the Pomeron-Pomeron — > x fusion vertex. 

Besides the pure exclusive process pp — > P + X+P: we will also study central diffractive inclusive x production 
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in which both of the incoming protons are destroyed, 

pp^X + x + Y. (2) 

Recall that for a heavy centrally produced state, such as a Higgs boson, the cross section for inclusive double- 
diffractive production is much larger than that for exclusive double-diffractive production OEIj 

find 3> CTcxcl- (3) 

The situation is quite different for the production of the relatively light x states. In this case we expect exclusive 
production to dominate. The reason is explained at the end of the next Section. Another interesting feature is 
that we can justify the use of perturbative QCD for inclusive double-diffractive x production by selecting events 
with large momentum transferred through the exchanged Pomerons (which can be measured as the tranverse 
energy flows, Eit, in the proton fragmentation regions). These Et^s provide the large scale which clearly allows 
the use of the perturbative QCD formalism. 



2 The Regge framework for double-diffractive x production 

Exclusive double-diffractive x production by Pomeron-Pomeron fusion is shown schematically in Fig. 2(a), where 
Pi,2± are the transverse momenta of the outgoing protons, 

^-^g%{pl^) V\p\^,pI^) sUpI^) (4) 



da 



(Ppi±d'^P2±dy 2^TT^ 

where (/at is the Pomeron-nucleon coupling and V is the Pomeron-Pomeron-x fusion vertex. The Pomeron 
trajectory is ap{t), with ti ~ —Pi±, and 

= ^ (5) 

where y is the cm. rapidity of the produced x meson, and 

M^^ = Ml + {pi^+P2i_f. (6) 

For simplicity we have neglected the real part of the Pomeron amplitude, assuming^ that Re/Im <C 1. The 
coupling is known from pp data, 

and ap{t) can be taken either from J/ip diffractive photoproduction data or from an analysis of pp data. The 
vertex V will be discussed later. 

For inclusive double-diffractive x production (Fig. 2(c)) we treat low- and high-mass proton dissociation 
differently. For low-mass dissociation, say TV N* , we just need to replace the coupling gN{t) by gN-^N' {t)- On 
the other hand the dissociation into a high-mass system is described in terms of the triple-Pomeron interaction 
with coupling g^p, see Fig. 2(c). The corresponding cross section is 



MfMida _ 1 2 , 2 ^TA2, 2 2 ^ . 2 ^ f Mf 



2 ]\j2 \ ap(0)-l 



dM^dMid^pi±d^P2±dy 2»TT 



Due to the condition 



Y9n{0) g3pipii_) V\pi^,pi^_) gspip: 



2±) 



X2S X /a;is\ 



(7) 



XIX2S = M^_L, (8) 

we see that the factor 

\ 2(qp(p? I )-l) 

gf) « .f.-H,;.„_ 



^In other words, we neglect the signature factors in Jl) and J7), assuming that ap{i) is rather close to 1. 
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Figure 2: (a) The Pomeron-Pomeron fusion diagram for double-difFractive exclusive production of a x(O^) 
meson, pp p + x+P] see (@J. (b) The rescattering or absorptive corrections to diagram (a), where the shaded 
region represents the soft pp rescattering corrections, leading to the suppression factor S'^. (c) The amplitude 
(squared) for inclusive double-diffractive x production; the final diagram shows how high-mass diffraction is 
given in terms of the triple-Pomeron interaction; see iQ. 
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and vice- versa. Thus the cross sections of both the exclusive and inclusive diffractive dissociation processes have 
the same energy and rapidity dependences, see JSJ. For fixed masses Mi^2 (and fixed Pi_l,P2_l and y) we have 

da cx s"p(p?i)+"^'(pL)~2, (10) 

However we have to include the absorptive corrections. The gap survival factor S'^ — that is the probability that 
the rapidity gaps will not be filled by secondaries produced in the soft rescattcring^ — is shown schematically 
by the graphs in Figs. 1(b) and 2(b). The corresponding amplitudes interfere destructively with the original 
bare amplitude. The elastic amplitude, shown by the shaded blob in Figs. 1(b) and 2(b), grows with energy and 
hence leads to a decrease of S^. The absorption is stronger for exclusive production (Fig. 1) since it is harder to 
transfer the loop momentum k± through the multiparticle chain formed by the diffractive excitations Mi and 

M2 [EI 3. 

First, we calculate the rapidity gap survival factors for exclusive and inclusive double diffractive Xc 
production using the formalism of Ref. ^Hl- If we were to neglect the slope of the Pomeron trajectory, ap — 0, 
then we find 

S'2(Tevatron) = 0.05 (0.15) (11) 

S'2(LHC) = 0.024 (0.10) (12) 

for exclusive (inclusive) production. When we include the Pomeron slope, ap = 0.1 GeV^ (which is consistent 
with the HERA J/ip photoproduction data ^HDi the above survival factors become 

S'^(Tevatron) 0.07 (0.16) (13) 

S'^(LHC) = 0.04 (0.11); (14) 
that is, in the Tevatron-LHC energy interval we may approximate the energy dependence by 

cx s-^ with S ^ 0.16 (0.1). (15) 

For inclusive production we select events with rapidity gaps A77 > 3: so ap acts on a smaller rapidity interval, 
and therefore is not so effective in enlarging the value of S*^ . 

Recall that the triple- Pomeron vertex g3p{0) is small numerically in comparison with the elastic Pomeron- 
nucleon coupling gAr(O); see, for example, ^Hl- On the other hand g^p is almost independent of the values of 
Pi± and P2±- Thus the suppression {gsp/gN)'^ is partly compensated by the larger interval of the transverse 
momenta pi± sampled in diffractive dissociation'*. To obtain a naive estimate of the ratio of the cross section 
of exclusive to inclusive x production, we first note that the ratio of cross section for single proton dissociation 
to that for elastic scattering satisfies |17[I13| 

R = ctsd/ctcI < 1/2. (16) 

Thus for the 'soft' process we expect 

CTincl/fTcxcl ^ < 1/4. (17) 

As a consequence, even without tagging the forward protons, we will mainly observe the exclusive pp ^ p + x+ P 
process simply by selecting events with rapidity gaps^. 

^See, for example, Ref. 1121 . and references therein. 

•^In impact parameter b± space, the elastic amplitude ylci(fex) arising from the inelastic channels, via unitarity, is given by the 
square of the inelastic matrix element, Aoi(6x) l^incl(^±)Pi and is therefore concentrated at smaller 6^ where the absorption is 
stronger; for more details see, for example, 1141 and references therein. 

''For the inclusive process, the momentum transfers pi± can be determined by experimentally measuring the transverse energy 
flows, EiT, in the proton fragmentation regions. 

^As we will discuss below, this is not true for heavy boson production, where we have large Sudakov T-factor suppression in 
order to ensure that the rapidity gaps survive against QCD radiation in the hard subprocess. Since the Sudakov suppression is 
much stronger for the exclusive amplitude for, say, Higgs production, we have crind ^ foxcl- 
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3 The perturbative QCD approach to double-difFractive x produc- 
tion 



To calculate the cross section for the central exclusive double-difFractive production of x(O^) states we use the 
formalism of Refs. [311213]. The amplitudes are described by the diagram shown in Fig.^a), where the hard 
subprocess gg x{^'^) is initiated by gluon-gluon fusion and where the second i-channel gluon is needed to 
screen the colour flow across the rapidity gap intervals. Ignoring, for the moment, the screening corrections of 
Fig. n^b), the Born amplitude is of the form [TH! 

T = A^^ [ ^.^/^t^l^^-p^ _ ^ f,{x,,x',,Ql^?■M)h{x2,x',,Ql,^^^■M), (18) 
Q\^(Q±-pi±Y{Q±+P2i.Y 

where the gg x(O^) subprocess is specified by jj] 

= SnTix gg)/M^, (19) 

and the polarisation factor^ 

^(x(0+)) = iQ±-Pi±) ■ (Q± +P2±). (20) 

Here T{x gg) is the width including the NLO corrections (that is the K factor). We assume the same 
NLO correction for the gg —>■ x vertex as for the X ~^ 99 width, which can be valid only within a certain 
approximation. 

The /g's in (|18|) are the skewed unintegrated gluon densities of the proton at the hard scale /i, taken typically 
to be Mx/2, with 

Qi = min I Q±,\{Q± - pi 
Q2 = mm I Q±,\{Q± + p2 
The longitudinal momentum fractions carried by the gluons satisfy 



(21) 



(22) 



Below, we assume factorization of the unintegrated distributions. 



U(x,x' ,0" ,iY-t) ^ jg(x,x' ,0" ,p')Fj,[i), (23) 

where we parameterize the form factor of the proton vertex by the form Fff(t) = cxp{bot) with bo = 2 GcV^^ 
|13|. In the domain specified by (|22|l the skewed unintegrated densities are given in terms of the conventional 
(integrated) densities g{x,Qf). To single log accuracy, we have [T^''. 

fg{x, x\ Qlfi^) = Rg-g^ {VT{Q^, m) xg{x, Q2)) , (24) 

where T is the usual Sudakov form factor which ensures that the active gluon does not emit additional real 
partons in the course of the evolution up to the scale, /i, of the hard process, so that the rapidity gaps survive. 
This Sudakov factor T is the result of resumming the virtual contributions in the DGLAP evolution. It is given 

by 



T(Q.,.)=exp -/ -^^l 



99 \ 

q 



dz . (25) 



Here, in analogy to |20| . we go beyond the coUinear approximation and in the T factor we resum, not just the 
single collinear logarithms, but the single soft log 1/(1 — 2;) terms as well. To a good approximation this can be 
achieved by taking the upper limit of the z integration in (|25(l to be 

A = . (26) 

kt + 0.62M^ ^ ^ 

^For x(0~) production the factor P(x(0~)) = ((Q± ^ Pi±) X {Q± + P2x)) ' "Oi where Hq is a unit vector in the beam direction; 
see, for example, 1181. 

^In the actual computations we use a more precise form as given by eq. (26) of Ref. Iljjl . 
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The square root in (|24|l arises because the bremsstrahlung survival probabiUty T is only relevant to hard gluons. 
Rg is the ratio of the skewed x' <C a; integrated distribution to the conventional diagonal density g{x,Q^). For 
a; <C 1 it is completely determined PT] , 

To compare expressions (0} and (the square of) ((TH)l . we recall, assuming Regge factorization, that the 
unintegrated gluon density can be written in the form 

fg{x,x',Q^fi^-t) = Rg a;^(Q''*) ^(Q^M^) ^^^^^^ (27) 

These factors x^, in the fg's in (|18|l . play the role of the Regge factors in (0J. Moreover the form 

factor FN{t) in (^5)) describes the t dependence of the coupling gwit), while the remainder of the integral in 
(O gives the product Vg%{0). 



4 The cross section for exclusive difFractive x production 

The uncertainties in the perturbative QCD predictions for the exclusive double-diffractive production of heavy 
(Higgs) states of mass 120 GeV were discussed in Ref. [20]. The uncertainties come from the infrared contribution 
to the Q± integral, (|18|l . from the choice of factorization scale, from the lack of precise knowledge of the gluon at 
low scales and small x, from deviations from formula H24|) for large-angle gluon emission and from the uncertainty 
in the calculation of the screening correction, S'^. For a Higgs 0+ state the total uncertainty was estimated to 
be a factor of almost 2.5, that is up to almost 2.5, and down to almost 1/2.5 the quoted value. In Ref. pil] 
it was also estimated that the uncertainties for producing a heavy 0^ state would be larger by almost another 
factor of 2.5. 

For the exclusive double-diffractive production of light x states the infrared uncertainties will be even larger. 
We computed the integral in 1)181) using the GRV94 parametrization H!]- the only available parton set which 
extends down to a very low scale (Q^ ~ 0.4 GeV^). To be safe, we choose a rather low energy, ^/s = 60 GeV, 
which corresponds to x ~ 0.05, where the GRV, MRST, CTEQ gluons are approximately equal; also where the 
pp total cross section is still approximately flat, which provides a simple normalization of the Pomeron-nucleon 
vertex. We find that half of the contribution to the amplitude (|18|) comes from the region Q± < 0.85 GeV. 
Clearly perturbative calculations are not justified in this domain. For low Q± the gluon propagator will be 
modified by non-perturbative dynamics; for example by the presence of GJ^j^G"^ condensates, effective gluon 
masses, confinement forces etc. It appears that these dynamics effectively suppress the low Q± contribution; 
for recent reviews see, for example, Ref. |^. Therefore the predictions below should be regarded as just an 
order-of-magnitude indication of the expected rates. If we evaluate amplitude (|18|l using GRV94HO partons 
[22 then we find, for y/s = 60 GeV, that the 'bare' cross section is 



da 



dpi±dp^±dy 



8 /ib/GeV'', (28) 



=P2±=0, y=0 



whereas if the contribution from the region Q± < 0.85 GeV is neglected then we obtain 2 /ib/GeV''. By the 
bare cross section we mean the prediction before including the rapidity gap survival factor . 

We emphasize that the diagram of Fig. 1(a) plays the dominant role when the transverse size of the Pomeron 
(that is the transverse separation of the x and x' gluons) is much larger than the size of the produced boson. 
This is the case for Higgs production^, but it is not so clear for Xc production. The other extreme is to assume 
that the transverse size of the Pomeron is relatively small, so that the Pomeron couples to each individual quark. 
This extreme corresponds to the 'additive quark model'. In this case the fusion vertex may be calculated from 
the Feyman diagram shown in Fig. 3. The coupling of the charm quarks to form the Xc meson is described in 
terms of the Xc meson wave function, 'tjj^r), for which we apply the standard non-relativistic formalism, see, for 
example, Refs. PSIOS] . 

Next we have to estimate the coupling gc of the charm quark to the Pomeron. For this, we can use the 
analysis of the data for charmonia production on nuclei [22] • The analysis [2H| showed that the cross section for 

* Indeed it has been shown 1241 that there is a negligible contribution from graphs where both gluons (forming the Pomeron) 
couple to the top-quark loop. 
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Figure 3: The diagram used to calculate the non-perturbative contribution to the exclusive process p+x+P- 
The amplitude, A^^, describing the Pomeron-Pomeron — > x subprocess, is controlled by the charm quark loop. 
It is evaluated using (j221)- 



J/^ on protons can be parametrised as^ 

= CTo (s/so)^, (29) 

with A = 0.21 and ctq = 3.6 mb, for sq = 100 GeV^. Thus the additive quark model estimate of the charm 
quark-proton cross section^" is cr(cp) ~ 2 mb, for s^p ~ 200 GeV^, which corresponds to -^/s = 60 GeV of 1^5)1 . 
Phenomenologically, the Pomeron-quark coupling may have a scalar or Dirac 7^ matrix form. The latter means 
that Pomeron exchange is similar to photon exchange. Thus it is reasonable to estimate the PP — x coupling 
using the known x ~* 77 decay width, r(x(0^) — > 77) ~ 3 keV |30| . However the analogy between the Pomeron 
and the photon cannot be exact. First, the Pomeron has even C-parity. Moreover the Pomeron is not a gauge 
field and so its 7^ vertex would not correspond to a conserved current |31| . Indeed, unlike x ~* 77i foi' PP ~^ X 
there are no graphs with a permutation of the two Pomeron vertices. The reason is that in Regge theory all 
interactions are strongly ordered in rapidity space. Thus we omit the second term in eq.(2) of Ref. |26| . 
If we follow the above procedure, then, for Xc(0"'"), we find that the PP — > x vertex, 1/ of l@J, is 

V ^ A^^ ^ gl ^ ±, (30) 
V6 rnc 



whereas the x ~* 77 decay amplitude is 



= \ (i'-'.^o) ^ (31) 



where gc is the Pomeron-c quark coupling, and a — ^JZjA.'Kmc 4' (0), where ij) (0) is the value of the first 
derivative of the bound-state wave function at the origin. In this way we can normalize the PP — > x vertex in 
terms of the known x ^ 77 decay width. We obtain 

WnaQED 

If we were to choose a scalar Pomeron-quark coupling, then we would obtain an App amplitude which is a 
factor 5/4 larger. However we prefer the 7^ form of the vertex, which conserves the s-channcl hclicity of the 
quark. 

If, now, we take cr(cp) = 2 mb, and normalize the A^^ amplitude as in (|32|l . then we find that the non- 

^The values of the cross section are in qualitative agreement with the analysis of Ref. 1291 . 

-'^'^Note, however, a similar analysis of 1/1 gives a somewhat larger value of (T(cp). 
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perturbative contribution -"^^ to exclusive double-difFr active Xc cross section is 



^^nonpert 



-1.1 fih/GeV^ X 5*2, (33) 



with ti ~ —Pi^- Note that the contribution of Fig. 3 is real (if we assume Re/Im^l), and therefore does 
not interfere with the imaginary amphtude of Fig.l. Thus we simply add these two contributions to the cross 
section. To be conservative, we only include the perturbative contribution of gluons with Q± > 0.85 GeV: that 
is 

^ - 2Mb/GeV%5^ (34) 

at y/s — 60 GeV; see the discussion below (j2Hl)- 

Now we need to consider the energy dependence. First we note that the data for J/^p photoproduction jl5| 
are described by a Pomeron with ap{0) — 1 ~ 0.2. The cc system in this process has a similar mass to the 
subprocess we are concerned with and so we take ap(0) — 1 = A = 0.2. The same value of A describes the 
behaviour of the JZ-^-proton cross section, and is also consistent with the small x behaviour found in the global 
parton analyses of deep inelastic scattering data. Adding together the above perturbative and non-perturbative 
contributions at y/s = 60 GeV, and using the appropriate energy dependence, s^'^ of lfTn|. we find the cross 
section for pp ^ p + Xc + P, at the Tevatron energy, ^/s = 2 TeV, is 



dy 



^ 130 nb. (35) 

y=0 



Here we have included the survival factor = 0.07 of the rapidity gaps, (|13|) . and integrated over the transverse 
momenta To obtain an approximate expression for the t dependence of the exclusive double-diffractive 

cross section, we neglect the possible difference in the form factors corresponding to Figs. 1 and 3, and the effects 
of the screening corrections on the pii_ distributions, which results in the behaviour F/v(ti)^F/v(t2)^ exp(Q! {ti + 

To obtain the total observable cross section at the Tevatron we integrate over an effective rapidity interval 
of Ay ~ 5 which gives 

CoxcKXc) = 650 nb. (36) 

This, more detailed estimate, is very close to our previous result 0. 

Note that the production of a (non-relativistic) 2++ quarkonium state may allow the discrimination between 
the two mechanisms shown in Fig. 1 and Fig. 3. As we already discussed (see Refs. j2| El for more details) in 
the on-shell approximation the Jz{gg) = amplitude vanishes. Then, the only contribution comes from the 
amplitude for jj^j = 2 quarkonium (2"^"'") production, arising from the violation of the = selection rule 
due to non-zero p±^ of the protons. Thus, we expect that the amplitude corresponding to Fig. 1 vanishes 
for forward 2^^^ quarkonium production. On the other hand there is no selection rule similar to Jz{gg) — 
for the general amplitude of Fig. 2 (though even P and C still hold). In particular, the amplitude of Fig. 3 
allows 2'^'^ production at any p± and, therefore, does not vanish as p± 0. Another difference between the 
contributions of Fig. 1 and Fig. 3, is that, in Fig. 3, a centrally-produced x(2++) meson has = 0, while for 
Fig. 1 non-relativistic 2++ quarkonium has \Jz\ = 2. 

For completeness we gather together all our predictions for the cross sections for double-diffractive x produc- 
tion in Table 1. The evaluation of the inclusive cross section is discussed in the next Section, and also in Section 
8. Comments on Xb may be found in Section 6. Although the cross sections are large, we note that the relevant x 
decay branching fractions are small. Traditionally, the experimentally favoured decay is Xc — *■ 7>//'0 ~^ 7/^^/^~ 
with branching fraction 6 x 10~^. On the other hand, Xc (tt+tt" + K^K~) has branching fraction 0.011, and 
this 2-prong configuration could prove more favourable. For Xb the corresponding branching ratios are not so 
well known, but it hard to imagine that the situation will be better. Recall that Xb decays dominantly into two 
gluon jets. Unfortunately this decay channel will suffer from a large background from QCD dijet production in 

^^We call it "non-porturbative" although, strictly speaking, the true underlying dynamics of the cr[cp) interaction is not well 
established. 
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70 
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Table 1: The cross sections (in nb) for exclusive and inclusive double-diffractive Xc and Xh production at the 
Tevatron and LHC. From the discussion in the text, we see that these predictions should be regarded as order 
of magnitude estimates. The accuracy is limited because of the low scale associated with x production. 

the same kinematical configuration. Using the results of Ref. 'S', and neglecting the K- factor in the background 
calculation, we estimate the signal-to-background ratio to be '--^ 5 x 10~'^/AAi^(GeV). Certainly we would 
expect that there is a sizeable (on the ~ 1% level) branching ratio for Xb(O^) 7 + T decay. But at the moment 
there is no experimental data on this channel. Therefore the situation regarding detection of the Xb(O^) in the 
central detector remains unclear. 

What are the uncertainties in our cross-section predictions? First, it is interesting to note that the results 
are almost independent on the choice of the slope, b = 26o + 2a Inl/x, of the pi±^ distributions. The reduction 
caused by a larger slope is largely compensated by an increased survival factor S"^ for the more peripheral 
reaction. For example, varying the slope 2&o from 4 to 6 GeV~^ changes the individual values of and iP' by 
approximately a factor of 2, but decreases j}? by only about 10-15% at the Tevatron and the LHC energies"'^^. 
On the other hand, the cross section would be 4 times smaller if we were to take a 'soft' Pomeron intercept, 
a_p(0) — 1 = A = 0.1, instead of 0.2, for the extrapolation from ^/s — 60 GeV to the Tevatron energy. Moreover 
the charm quark cross section may be larger (see the footnote below (|29|l 'l. and there may be an additional 
contribution from the low domain, Qj^ < 0.85 GeV, of the gluon-induced amplitude lfTH|l . Thus the cross 
section estimate of H36|l could be a factor of 4 higher. 



5 The rates for inclusive double-diffractive x production 

The cross section for double-diffractive inclusive production is calculated is an analogous way to the exclusive 
one. Indeed, it may be considered as 'exclusive' x meson production in parton-parton collisions. The flux 
of initial partons is given by the known parton distributions inside a proton, while the unintegrated skewed 
distributions (inside an individual) parton, /g, are calculated in terms of the non- forward BFKL amplitudes. It 
is important to note that here there is no form factor for point-like parton scattering and that there are sizable 
contributions over a much larger range of transverse momenta p^^, with i — 1,2. Therefore we keep just the 
leading-order contributions; that is, the integrals over Qj^, and take logarithmic forms. In particular 

'^^^-^ takes the form ^ (37) 



(38) 



Qi(Q± -Pi±)2(Q± +P2i)2 Q\ pIi_ 

so that the final integrals in the expression for the cross section become 

(PQi_ (PQ'i^ (Ppi± (Pp2± 
Ql Q'l P\^ pI^ ' 
see Refs. [321 Ej for details. 

As mentioned before, for a heavy centrally produced state, such as a Higgs boson, the cross section for 
inclusive double diffractive production is predicted to be much larger than that for exclusive production. The 
reasons are as follows: 



^^For the same reason the cross section predictions for exclusive double-diffractive Higgs production are practically independent 
of the value taken for the slope h. 
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(a) due to the absence of proton form factors in the inclusive process there are sizeable contributions to the 
dp\lp\ integral with Q_l < p_L < Mh/'2., 

(b) the T-factor suppression is weaker for inclusive production, as now T{p±, Mh) > T{Q±, Mh), 

(c) for inclusive production there is no P-even ~ selection rule [331 1^ 12 El > which suppresses^^, for 
example, 0^,2+ (non-relativistic) quarkonium exclusive production pi 151 ITl]. 

On the other hand for the relatively light x states the typical values of p± transferred through the Pomeron, 
and Q± of the internal gluon (inside the Pomeron) are of the same order. As a consequence inclusive production 
is not strongly enhanced in comparison with exclusive production. Insight into this result can be obtained from 
the discussion leading to (|17|l . It is therefore not surprising that, using the perturbative QCD approach of the 
previous section, we estimate that 

ainci(Xc(0+)) ^ 70 nb (39) 

at the Tevatron energy y/s = 2 TeV, which is about 10 times less than the exclusive double-diffractive cross 
section CTcxcI — 650 nb. It is interesting to note that the perturbative QCD estimate reflects the smallness of 
the triple-Pomeron vertex. The low QCD estimate arises from the small gluon density (at low x and small 
scale) obtained in the global parton analyses of deep-inelastic and related data. As we have already mentioned, 
recently the CDF collaboration 10 have reported their first results on J/ip^^ central production, with rapidity 
gaps on either side, and with a mass consistent with that of Xc- These results are still preliminary, but the 
observed cross section is compatible with the expectations of Ref. assuming that the background issues are 
well understood and turn out to be small. With more data to come, and with a study of the background, it 
will be possible to obtain more detailed information on this process. 

In going from the Tevatron to the LHC energy, y/s — 14 TeV, the double-diffractive inclusive cross section 
(for a fixed rapidity gap interval), Uinci, is expected to grow as s^, with A ~ 0.2. Now we have to integrate 
over the rapidity gap intervals A77 > 3, that is over the masses dMf/Ml and dM2/M2. In the limit of small 
fixed masses Mi, M2 the integrand grows faster - as s^^ - similar to the elastic cross section. On the other 
hand the rapidity gap survival factor S*^ decreases as s~°-^ ^31, see (|15|l . To be explicit, we find [T3| 

5*2 (LHC)/^^ (Tevatron) ~ 2/3. 

We see that the energy dependence of the cross section and the behaviour of the survival factor partially 
compensate each other, so 

^^(LHC) ~ 2.5 (Tevatron), (40) 

dy dy 

for 2/ = 0. At the LHC, a larger rapidity interval is available. Ay ~ 8, which leads to an increased cross section 

(Tinci(Xc(0+)) ^ 200 nb. (41) 

For the exclusive process, the cross section da/dy is proportional to s^^. In fact the QCD estimate gives 
(Tcxci(Xc(0''')) — 3000 nb at the LHC. These predictions are collected together in Table 1. 



6 Predictions for double-diffractive Xb production 

The perturbative QCD predictions of the cross sections of double-diffractive Xfc production are more justified 
than those for Xc, on account of the larger Xb mass. For this reason we only evaluate the perturbative amplitude 
l|18|l . which corresponds to the diagrams of Fig.l. Moreover, as before, we keep just the contribution from the 
Q± > 0.85 GeV integration region. This is found to give about | of the total amplitude H18I) obtained using 
GRV94HO partons 122. The energy dependence of the bare cross section is driven by the small x behaviour 
of the gluon. However we use again A — 0.2, which is much more consistent with the recent partons sets (for 
example EEl) than the GRV94 set. 

Another uncertainty is the absence of an experimental measurement of the gluonic width, T{xb{0^) gg)- 
Potential models |371 predict a bare width of ro(x6 — ^ gg) = 1.3 — 1.9 MeV, which should be multiplied by the 

^^The exclusive double-diffractive production of vector and axial vector states is suppressed as a consequence of the Landau- Yang 
theorem 1351 for massless gluons. 
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NLO correction factor K = (l+Q.Sas/Tr) [3HI- This is some 5 times larger than the prediction ro(xfc — > gg) ~ 350 
keV obtained from QCD lattice calculations IMI ■ For the predictions of the cross sections of double-diffractive 
Xb production at the Tevatron and LHC, which we show in Table 1, we have taken r{xb gg) — 1-3K MeV. 
Note that these predictions are larger than those given in Rcf. where the lattice estimate of the width was 
used^'*; see also Ref. [5]. 

7 Implications of the larger exclusive rate for Xc production 

The QCD predictions of the double-diffractive exclusive and inclusive Xc production show that the former 
process dominates. Indeed at the LHC we estimate 

Crcxcl/o-incl ^ 10. (42) 

It means that by just selecting events with two rapidity gaps we will observe mainly the exclusive process, even 
without requiring the forward protons to be tagged! 

One topical application would be to search for the exclusive production of the new X(3872) charmonium 
state observed recently by the Belle and CDFII collaborations 42 . If this state were a radial excitation of the 
Xc(O^) meson, then it would be observed with a cross section 

a{pp -^p + X+p) - 100(500) nb (43) 

at the Tevatron (LHC). Of course, looking for a specific decay X V^tt+tt^ /i+/^^7r+7r^ will reduce the 
signal by about 1000, but still the cross section is large. 

Since the X production rate is proportional to the X ^ gg decay width, a measurement of the pp p+X+p 
cross section can give an estimate of the X ^ gg coupling. To reduce the (large) theoretical uncertainties it 
would be better to compare the rates for X(3872) and Xc(3415), as the gg width is known for the latter state^^. 

If X(3872) were a £'D*-molecule (see, for instance, 03]), then the exclusive signal would probably not be 
seen^^. Also note that, due to the P-even, = selection rule^'', the exclusive cross section is suppressed 
for the double-diffractive production of 0^, 2+ or C-odd bosons. Thus by observing the X(3872) state in the 
exclusive process we would conclude that it is a C-even, and most probably a 0+, particle. 



8 Parity determination 

With sufficient luminosity we can use the inclusive process, to study 0^ or 2++ production by Ponieron- 
Pomeron fusion. If the transverse momentum p± transferred through the Pomeron is small, p± <C Q±, then 
the inclusive process is controlled by the same P-even, Jz = selection rule as exclusive production. However 
if we select events with large pi± then the production of P-odd and Jz — ±2 states become possible. Note 
that large p±^s can be measured as the transverse energy flows, pix = — and p2i_ — E2_l, in the proton 
fragmentation regions. Moreover the major part of the energy flow is carried by one jet (with lowest rapidity 
|y| in the centre-of-mass frame). Indeed at lowest order this jet carries the entire Ei±. These measurements 
allow a study of the dependence, where (f> is the azimuthal angle between the energy flows Ei± and E2±. The 
(p distribution depends on the parity of the produced system |20j : 

da/d(f> oc 1 -I- cos2(/i for a natural parity state, P = (—1)'' (44) 

^*Note, however, that the calculations of the matrix elements in Ref. 1391 were based on the lattice NRQCD results of Ref. I4UI . 
which used quenched approximation. As was shown in f41| . this approximation underestimates the NRQCD matrix elements. Thus 
the values of the widths, in the case of three light quark flavours, should be higher than the quenched results. 

^^The effective gg^^ luminosity, C{M^), is essentially flat in this mass interval. 

^^In this case there should be some contribution from the D-meson loop, analogous to the charm quark loop in Fig. 3, but it is 
expected to be small. 

^'^Unfortunately the Jz = selection rule is less precise at these low scales. The admixture of |Jz| = 2 states may be up to 
20-30%. However, it is encouraging that expectations based on this rule appear to be in good qualitative agreement with the 
available data on double-diffractive meson production; sec 14 for details. 
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E^>3 




E±>7 












500 


20 






35 


5 




LHC 










6000 


400 


40 




500 


100 


30 



Table 2: The inclusive double-diffractive cross sections (in pb) for various choices of the cut on the transverse 
energy flows, Ei± (in GcV). To account for the large gluon virtuality, which occurs for large Ei±, we follow 
a procedure similar to that used in Ref. for the calculation of the amplitude for the decay of x(O^) iiito 
virtual photons. 

whereas 

da/dcj) cx 1 — cos2(/) for an unnatural parity state, P = —(—1)'^. (45) 

We emphasize, that since we consider rather large E±, these angular distributions are less sensitive to soft 
rescattering. 

To get some idea of the event rates, we present the cross sections for pp ^ X + x(0+) + F in Table 2, for 
different choices of the E± cut. If we take Ei± > 3 GeV at the Tevatron, then a'inci(Xc(0^)) ^ 0.5 nb, whereas 
for Ei± > 7 GeV at the LHC we estimate crinci(Xc(0^)) ^ 40 pb. Note that the reliability of the predictions 
with large E± cuts is better, since we do not enter the infrared domain. These are rather large cross sections, 
but recall that the relevant decay branching fractions are small, see Section 4. 

Finally, note that the ratio of the 0++ and 2++ production cross sections is 

^inci(2++) ^ 5 r(2++) 

ainci(0++) 4 r(0++)' ^ ' 

where we have neglected the mass difference between the M(2++) and Af(0++) mesons. 



9 Conclusions 

We find that both the Regge formalism and perturbative QCD predict essentially the same qualitative behaviour 
for the central double-diffractive production of 'heavy' Xc(0^^) and Xb{0^'^) mesons. Due to the low scale, 
there is a relatively small contribution coming from the process in which the incoming protons dissociate. 
Therefore simply selecting events with a rapidity gap on either side of the x, almost ensures that they will come 
from the exclusive reaction, pp p + x ^ P- 

We evaluated the expected double-diffractive cross sections and demonstrated that they are sufficiently large 
for both Xc and Xb meson production to be observed. Since the rapidity-gap survival factor, S*^, decreases with 
energy, the cross section da/dy at the LHC energy only exceeds that at the Tevatron energy by a factor of 
about 3. Our cross section predictions should be regarded as only order of magnitude estimates, because of 
the low scale, nevertheless these processes can be very informative. For example, the observation of the new 
charmonium state X(3872) in the exclusive process pp p + X + p would be a strong argument in favour 
of its quantum numbers being J^^ = . Moreover, since the exclusive cross section is proportional to the 
gluonic width, it will be possible to measure the width T{X gg) by comparing the rates of exclusive Xc and 
X production. While there is a sizeable uncertainty in the predictions for the overall rates of double difffractive 
X production, we stress that the kinematic distributions, for example the x transverse momentum distributions, 
should be more reliable. 

Although exclusive x production is expected to dominate, the event rates should be large enough to select 
double-diffractive dissociative events with large transverse energy flows in the proton fragmentation regions. 
Such events are particularly interesting. First, in this case, the large value of E±^ provides the scale to justify 
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the validity, and the reasonable accuracy, of the perturbative QCD calculation of the cross section. Next, by 
measuring the azimuthal distribution between the two E± flows we can determine the parity of the centrally 
produced system. 

An interesting extension of the exclusive double-diffractive approach, would be to observe central open bb 
production; namely 6, b jets with p± <^ mf,. Again, this would put the application of perturbative QCD on a 
sounder footing. It would allow a check of the perturbative formalism, as well as a study of the dynamics of bb 
production. 
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